Shadow on the wall cast by an Abrikosov vortex 
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At the surface of a d-wave superconductor, a zero-energy peak in the quasiparticle spectrum 
can be observed. This peak appears due to Andreev bound states and is maximal if the nodal 
direction of the d-wave pairing potential is perpendicular to the boundary. We examine the effect 
of a single Abrikosov vortex in front of a reflecting boundary on the zero-energy density of states. 
We can clearly see a splitting of the low-energy peak and therefore a suppression of the zero-energy 
density of states in a shadow-like region extending from the vortex to the boundary. This effect is 
stable for different models of the single Abrikosov vortex, for different mean free paths and also for 
different distances between the vortex center and the boundary. This observation promises to have 
also a substantial influence on the differential conductance and the tunneling characteristics for low 
excitation energies. 

PACS numbers: 74.45.+c,74.20.Rp,74.25.-q 



Today, there is common agreement that most high- 
T c superconductors exhibit <i-wave symmetry. An im- 
portant characteristic of d-wave superconductors is the 
possible existence of Andreev bound states at its surface 
IlLlafcl- This increase of the local zero-energy quasiparti- 
cle density of states at the surface can clearly be observed 
in the differential tunneling conductance as a pronounced 
zero-bias conductance peak ("H- [3- (Z3| ■ For specular bound- 
aries, this peak reaches maximum height, if the surface is 
perpendicular to the nodal direction of the ei-wave. The 
effect shrinks, if the orientation is changed 0, 01- For 
an angle of 45 degrees between the nodal direction and 
the surface the bound states have vanished completely. 
However, it has been pointed out, that for rough surfaces 
a similar shape of the zero-bias conductance peak is ob- 
tained which is independent of the boundary orientation 
0- If a magnetic field is applied, the spectral weight of 
the zero-bias peak decreases and a splitting of the peak 
is observed HUGH- 

In this Letter, we study a single Abrikosov vortex in 
front of a specular surface, and we investigate the effect 
of the vortex on the local quasiparticle density of states 
along the boundary. All interesting phenomena of this 
problem are described within the quasiclassical theory, 
which is valid if the coherence length is much larger than 
the Fermi wavelength. To calculate the local density of 
states in the vicinity of the boundary it is necessary to 
find numerically stable solutions of the Eilenberger equa- 
tions 0, ^3 that fulfill the appropriate boundary condi- 
tions at the specular surface. For this purpose we use the 
Riccati parametrization of the quasiclassical propagator 
fbil. Along a trajectory of the kind r(x) = r + x v F the 
Eilenberger equations of superconductivity reduce to a 
set of two decoupled differential equations of the Riccati 
type for the functions a(x) and b(x) 

hv F d x a(x) + [2e„ + A t a(a;)] a(x) - A = 
hv F d x b(x) - [2e„ + Ab(x)] b(x) + A 1 = (1) 



where ie n = ie n +v F ■ -A are shifted Matsubara frequen- 
cies. For the simple case of a cylindrical Fermi surface 
the Fermi velocity can be written as 



vp = Wf(ei cos 9 + e2 sin( 



(2) 



The 9- and r-dependence of the pairing potential A can 
be factorized in the form 



A(r, 



(3) 



For a d x 2_ y 2 -wave superconductor the symmetry func- 
tion takes the form = cos(2(9) and for the s-wave 
symmetry it becomes constant = 1. For e n > the 
Riccati equations have to be solved using the bulk values 
as initial values at x = ±00 



a(— 00) = 
b(+oo) = 



A(-oo) 
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The numerical solution of the Riccati equations can be 
done with minor effort and leads to rapidly converging 
results. For the calculation of the local density of states 
the imaginary part of the quasiclassical propagator has to 
be integrated over the angle 9 that defines the direction 
of the Fermi velocity. In terms of a and b we have 
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where E denotes the quasiparticle energy with respect 
to the Fermi level and S is an effective scattering pa- 
rameter that corresponds to an inverse mean free path. 
As it is well known, the localized zero-energy state in a 
d a; 2_ 1 ,2-superconductor at a specular 110-boundary arises 
from the sign change in the pairing potential on a quasi- 
particle trajectory that is reflected at the surface. The 
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outgoing particles are Andreev reflected at this potential 
step and interfere with the incident quasiparticles. This 
interference process leads to stable zero-energy trapped 
states in the vicinity of the boundary, called Andreev 
bound states. The same sign change in the order para- 
meter is found on a trajectory that passes near the cen- 
ter of a vortex and therefore leads to similar localized 
Andreev bound states inside the vortex core. The sup- 
pression of the amplitude of the pairing potential around 
the vortex center gives only a small quantitative correc- 
tion in the calculation of the trapped state as we already 
pointed out before Q| • The influence of the boundary for 
anisotropic superconductors is included within the quasi- 
classical theory if the nonlinear boundary conditions for 
the quasiparticle propagator are applied [lallal- For the 
Riccati parametrization a substantial simplification oc- 
curs and an explicit solution of the nonlinear boundary 
conditions can be found [l7j . 

In the following we assume a totally reflecting surface 
where the transparency T equals zero while the reflec- 
tivity 1Z becomes one. In this special case the boundary 
conditions reduce to 

Oout = «in and °in = 6 out (6) 

In the next step we try to find an appropriate model 
that describes the pairing potential associated with the 
single vortex in front of the reflecting boundary. With 
the condition that there are no currents flowing across the 
boundary we have to find a phase of the order parameter 
where the phase gradient vanishes perpendicular to the 
boundary. In analogy to the classical boundary value 
problem of electrostatics with a point charge in front of 
a conducting surface we introduce an image vortex on 
the opposite site of the reflecting boundary. The pairing 
potential around a vortex at position ry can be written 
as (see also [T^ l 

*(r) = /(|r-rv|)e l *M (7) 

The function f(\r — ry|) characterizes the amplitude of 
the pairing potential of the single vortex. Since we con- 
sider a vortex-antivortex pair, the phase $(r) is given as 

$(r) = arg(r — ry) — arg(r — fy) (8) 

The location of the image vortex behind the boundary is 
defined by fy — ry — 2n(n, ry). The normal unit vector 
is given as h = l/^/2(e 1 + e 2 ) for a 110-boundary. 

In the following, the origin O of our coordinate sys- 
tem is placed at the boundary right between the vor- 
tex and the image vortex. The y-axis is orientated par- 
allel to the boundary, xy denotes the vortex position 
on the x-axis and also measures the vortex to boundary 
distance. Furthermore it is useful to introduce the co- 
herence length £ = as the unit of a general length 



scale. We performed calculations of the local density 
of states for both a model pairing potential modulus 
f(\r — ry|) — tanh(|r — ry|/£) and a constant modulus 
/(|r — TV |) = 1. The latter corresponds to a pure phase 
vortex. The results according to both models show only 
small quantitative differences. In particular, the main 
qualitative effect we want to present here, the shadow on 
the zero-energy density of states, exists independently. 
Thus, we will restrict our following considerations to the 
simpler second model. 

In Fig.^we show the zero-energy local density of states 
in the vicinity of the reflecting boundary. In the up- 
per part of the image, the local density of states is dis- 
played as a three-dimensional surface, in the projection 
below we show the same quantity as a density plot. A 
phase vortex is situated at a distance of two coherence 
lengths xy = 2£ from the surface. We assume a dis- 
symmetry of the order parameter and set the boundary 
with an angle of 45 degrees to the &-axis of the crystal. 
First we notice, that the fourfold symmetry of the local 
density of states around the vortex is broken. Then we 
also observe a shadow-like suppression of the zero-energy 
density of states in a triangular region between the vortex 
and the boundary. 

The picture for an s-wave superconductor is totally 
different. As shown in Fig. the vortex has only little 
influence on the boundary density of states. Here, the 
high zero-energy density of states of the vortex rather il- 
luminates the boundary. The small shadow on the right 
hand side of the vortex and the slight line between the 
vortex and the boundary is due to quasiparticles with 
an inclination angle of 90 degrees, that are trapped be- 
tween the reflecting boundary and the potential step in 
the vortex core. Below, we will focus on a d x 2_ y 2-w&ve 
superconductor with an 110-boundary. More details of 
the s-wave superconductor and a e?-wave superconduc- 
tor with an arbitrary orientation of the boundary will be 
discussed in a following work. 

In Fig. 13 we show the zero-energy density of states 
along the 110-boundary. The different curves correspond 
to different vortex positions xy. The calculations have 
been done using a value of 8 = O.lAo. For increasing 
vortex to boundary distances xy we find a decrease of 
the shadow depth, while the width increases visibly. Ap- 
parently, the shadow effect exists for a wide range of vor- 
tex to boundary distances xy even larger than 10£. For 
distances smaller than one coherence length £ a selfcon- 
sistcnt calculation of the pairing potential around the 
vortex might become necessary. 

In the following we want to obtain a more detailed 
impression of the vortex shadow at the point O. We 
also want to study the influence of the effective quasi- 
particle scattering parameter 5 introduced in Eq. (pj on 
the local density of states. Without vortex, the value of 
the zero-energy density of states at the boundary shows 
a sensitive dependence on the decoherence parameter d. 
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FIG. 1: Zero-energy local density of states of a d a ,2_ H 2-wave 
superconductor. A reflecting 110-boundary is situated at the 
back of the image, the single Abrikosov vortex can be found 
at a distance of two coherence lengths xy = 2£ in front of 
it. The Andreev bound states at the boundary are seen as 
the bright regions along the y-axis, and as the high values in 
the three-dimensional plot, respectively. The localized states 
in the vortex center are seen as the peak in the upper image 
and as the bright spot in the image below. One can clearly 
identify the distinct shadow that emanates from the vortex 
center to the boundary. An effective scattering parameter of 
S = 0.1 Aq has been used for the calculations. 



With decreasing 5 the value of the zero-energy peak at 
the boundary increases rapidly. In Fig. 0] we show the 
zero-energy density of states at the point O as a function 
of the vortex to boundary distance xy for different values 
of 6. The curves are normalized to the particular values 
of the zero-energy boundary density of states without a 
vortex or far away from the vortex center. We find that 
both the range and the relative depth of the shadow in- 
crease, if we decrease the scattering rate S. We want 
to point out, that this shadow effect can be observed in 
clean superconductors with a long mean free path as well 
as in superconductors with higher scattering rates. 

In order to explain the suppression of the local zero- 
energy density of states at the surface, we now concen- 
trate on a given point in the shadow region. In order to 
find the quasiparticle spectrum there, the angular inte- 
gration in Eq. JSJ has to be done. For each angle, the in- 




FIG. 2: Zero-energy local density of states of a s-wave super- 
conductor. Again the reflecting boundary is situated at the 
back of the image, while the single Abrikosov vortex can be 
found at a distance of one coherence length xy — £ in front of 
it. The high zero-energy density of states of the s-wave vor- 
tex slightly illuminates the boundary. An effective scattering 
parameter of 8 — 0.1 Aq has been used for the calculations. 



tegrand corresponds to the contribution of a quasiparticle 
trajectory with the direction specified by 6. Due to the 
phase gradient of the order parameter, the energy "seen" 
by a quasiparticle flying along a trajectory is shifted. Ad- 
ditionally, this shift itself changes locally along the trajec- 
tory. Thus, for most of the angles, the Riccati-equations 
(Eq. QJ) are not evaluated at zero-energy. This is suf- 
ficient, however, to miss the sharp zero-energy peak of 
the bound state at the surface. As a consequence, the 
zero-energy density of states is reduced. In the quasipar- 
ticle spectrum the spectral weight of the bound states 
is shifted from the Fermi level towards higher energies. 
This effect is similar to the splitting of the zero-bias peak 
due to surface currents 0. 

In Fig. [S] we show the local density of states at the 
point O for different vortex to boundary distances xy as 
a function of energy. If the Abrikosov vortex is placed in 
the vicinity of the boundary we observe a distinct split- 
ting of the zero-energy peak. With increasing distance 
between vortex and boundary the splitting is reduced. At 
xy = 10£ a splitting is no longer visible, while the height 
of the zero-energy peak is still considerably reduced. 

The strong reduction of the zero-energy density of 
states at the 110-boundary of a d x 2 _ y i -superconductor 
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FIG. 3: Zero-energy density of states of a d x 2_ y 2-wecve super- 
conductor along a 110-boundary for different vortex positions 
xv- The curves are normalized to the local density of states 
at the boundary without vortex (xv — > oo). The effective 
scattering parameter is chosen to be S = O.IAq. 



FIG. 5: Quasiparticle spectrum of a d x 2_ y 2-wa,ve supercon- 
ductor taken at a 110-boundary at the point adjacent to the 
vortex center (C). The curves are calculated for different vor- 
tex distances xv from the boundary. The effective scattering 
parameter is chosen to be S = O.IAq. 
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FIG. 4: Zero-energy quasiparticle density of states of a 
d x 2_ y 2-wa.ve superconductor at the boundary as a function 
of the vortex position xy. The different curves correspond 
to different mean free paths related to the inverse of 8. The 
DOS of each curve is normalized on the corresponding bound- 
ary density of states without vortex. 



has of course an important influence on the zero-bias 
anomaly in the tunneling conductance. Even in zero 
magnetic field vortices can remain in the high T c - 
materials by pinning defects. In the vicinity of grain 
boundaries these pinned vortices will play an important 
role for the grain boundary tunneling due to the reduc- 
tion of the zero-energy density of states. In a more de- 
tailed work we will also discuss the influence of a single 
vortex on the local density of states in the vicinity of 
a rough surface and we will consider several interesting 
boundary geometries apart from the flat surface. 

S. G. is supported by the 'Graduicrtcnforderungs- 
programm des Landes Baden- Wiirttemberg'. C. I. is 
grateful to the German National Academic Foundation. 
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